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1.0   INTRODUCTION 

Infinite  planes  of  gamma-ray  sources  are  of  interest  in  fallout 
radiation  studies.   The  energy  distribution  of  the  gamma-ray  flux  arisii 
from  such  a  source  is  usually  measured  using  a  scintillation  spectrometer 
systcm  consisting  of  a  light-producing  crystal,  photomultiplier  (I'M)  tube, 
associated  amplifiers  and  power  supplies,  and  a  multichannel  pulse-height 
analyzer.   The  accumulated  pulse-height  spectrum  represents  the  response 
of  the  system  to  the  gamma  rays  striking  the  crystal.   The  primary  infor- 
mation desired  from  such  a  spectrum  is  the  energy  distribution  of  the 
gamma-ray  flux  at  the  crystal.   In  order  for  such  an  analysis   to  be  pos- 
sible, the  response  of  the  spectrometer  system  to  photons  of  known  energy 
and  intensity  must  be  measured. 

The  energy  distribution  of  the  incident  flux  may  vary  with  direction. 
Consequently,  measurements  are  usually  made  with  a  collimated  detector, 
which  limits  the  measurements  to  a  specified  direction.   It  should  be  noted, 
however,  that  the  collimator  itself  influences  the  system  response. 

Differences  between  the  true  energy  distribution  and  the  pulse-height 
distribution  are  caused  by  two  types  of  distortion,  detection  distortion 
and  electronic  distortion.   The  distortion  in  the  detection  is  caused  by 
the  statistical  nature  of  the  gamma-ray  interactions  in  the  crystal  and  the 
scintillator  variance.   The  electronic  distortion  is  a  result  of  the  var- 
ious staged  processes  required  to  convert  the  light  output  of  the  crystal 
into  a  voltage  pulse,  and  the  subsequent  transfer  of  the  voltage  pulse  from 
the  output  of  the  PM  tube  to  the  multichannel  analyzer  (MCA) .   If  the  re- 
sponse of  the  system  to  monoenergetic  photons  is  measured,  one  has  in  effect 
measured  the  distortions  associated  with  photons  of  that  particular  energy. 
Once  these  measurements  are  completed  for  several  different  photon  energies, 


the  analysis  of  a  complex  experimental  spectrum  reduces  to  a  mathematical 
problem  of  determining  the  set  of  energies  and  intensities  which  best 
represent  the  experimental  spectrum  (3,  15,  12,  4). 

Typical  problems  encountered  are  the  measurement  of  the  angular 
distribution  of  barrier-transmitted  radiation,  or  the  angular  distribution 
of  gamma-ray  photons  reflected  from  a  barrier.   In  order  to  determine  the 
energy  distribution  of  the  flux  measured  in  these  types  of  experiments,  it 
is  necessary  to  know  the  system  response  to  a  uniform  flux  of  known  energy 
distribution.   The  basic  requirement  is  that  the  flux  be  constant  across 
the  collimator  aperture.   In  order  to  produce  a   flux  which  satisfies  this 
requirement,  an  infinite  plane  isotropic  source  is  simulated.   The  simula- 
tion may  be  accomplished  by  a  point  source  superposition   technique.    The 
use  of  the  KSU  collimator-spectrometer  system  to  measure  energy  and  angular 
distributions,  which  are  to  be  analyzed  using  the  response  functions  pre- 
sented in  this  work,  imposes  no  geometrical  restrictions.   The  only  restric- 
tion is  that  the  photon  flux  must  be  uniform  across  the  collimator  aperture, 

This  study  is  concerned  with  experimentally  measuring  the  response  of 
the  Kansas  State  University  spectrometer  system  to  photons  from  nine  mono- 
energetic  gamma-ray  emitters,  i.e.,  Pb-210  (47  keV) ,  Cd-109  (88  keV) , 
Ce-139  (166  keV) ,  Hg-203  (279  keV) ,  In-Sn-113  (392  keV) ,  Sr-85  (513  keV) , 
Cs-Ba-137  (662  keV) ,  Mn-54  (840  keV) ,  and  Zn-65  (1114  keV) . 


2.0  THEORY 

2.1  Detection 

The  response  of  the  spectrometer  system  is  primarily  a  function  of 
the  detection  process.   The  system  under  consideration  employs  a  Nal(Tl) 
scintillator  crystal  connected  to  a  ten  stage  PM  tube.   The  height  of  the 
voltage  pulse  at  the  output  of  the  PM  tube  is  proportional  to  the  amount 
of  energy  deposited  in  the  crystal  detector,  (3,17,4).   This  proportional 
nature  allows  the  use  of  a  multichannel  analyzer  which  stores  the  voltage 
pulses  in  a  magnetic  memory  unit  according  to  their  height;  thus   the 
accumulation  of  a  pulse-height  spectrum  is  realized.   The  observed  dis- 
tribution represents  the  total  response  of  the  system,  as  arranged,  to  the 
gamma  radiation  being  studied.   Figure  1  shows  a  typical  response  curve 
for  a  monoenergetic  gamma-ray  source  with  characteristic  X  ray,  indicating 
the  main  features  of  the  distribution.   The  shape  of  the  spectrum  is  a 
consequence  of  the  photon  interactions  in  the  crystal.   The  source  photons 
will  interact  with  matter  in  three  principle  ways,  photoelectric  absorp- 
tion, Compton  scattering,  or  pair  production. 

In  photoelectric  absorption  the  entire  energy  of  the  photon  is 
transferred  to  an  electron  in  one  collision.   The  Compton  effect  is  char- 
acterized by  the  fact  that  only  a  portion  of  the  photon  energy  is  trans- 
ferred to  an  electron  in  a  collision.   The  amount  of  energy  lost  by  the 
photon  is  a  function  of  the  scattering  angle,  as  shown  by  equation  (2.1), 
where  E  and  E1  are  respectively  the  initial  and  final  photon  energies. 
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The  photon  exit  angle,  <J>,  is  measured  with  respect  to  the  incident 

direction;  m  is  the  electron  rest  mass, 
o 

The  largest  possible  energy  loss  occurs  when  the  photon  undergoes 
a  head-on  collision  with  an  electron  and  scatters  180  .   The  smallest 
energy  loss  occurs  when  (J>  =  0  ,  i.e.,  E1  =  E.   The  Compton  distribution 
begins  at  zero  energy,  since  the  recoil  electron  has  zero  energy  for 
<J>  =  0  .   The  maximum  pulse  height  of  the  Compton  distribution  occurs  when 
$  =  180  ,  where  the  recoil  electron  has  its  maximum  possible  energy.   This 
maximum  pulse  height  is  referred  to  as  the  Compton  edge. 

Pair  production  is  a  threshold  type  of  interaction  between  the  inci- 
dent gamma  ray  and  the  nuclear  forces  surrounding  the  nucleus.   A  gamma 
ray  may  annihilate  itself  in  the  nuclear  force  field  and  produce  a  positron- 
electron  pair.   The  rest  mass  of  each  particle  is  0.511  MeV,  thus  the  pho- 
ton must  have  an  energy  of  at  least  1.022  MeV  in  order  to  produce  the  pair. 
Any  excess  energy  is  given  to  the  pair  as  kinetic  energy.   The  electron 
and  the  positron  do  not  travel  far  from  the  production  point  before  the 
positron  recombines  with  an  electron  and  annihilates.   The  annihilation 
results  in  the  production  of  two  0.511  MeV  photons  traveling  in  exactly 
opposite  directions.   Either  or  both  of  these  photons  may  interact  in  the 
crystal  by  Compton  or  photo  absorption  mechanisms. 

In  general   it  may  be  said  that  the  lower  energy  region,  beloi'/  about 
300  keV,  is  dominated  by  the  photoelectric  effect.   The  mid-energy  region 
between  300  keV  and   4  MeV  is  dominated  by  the  Compton  effect.   At  energies 
much  above   4  MeV  pair  production  becomes  the  dominant  effect.   These 
trends  are  illustrated  in  Figure  2,  which  shows  the  photon  interaction 
cross  sections  as  a  function  of  incident  energy,  up  to  about  2  MeV. 

Nal(Tl)  is  a  luminescent  material.   Gamna  rays  entering  the  crystal 
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detector  cause  ionization  and  excitation  of  the  Nal  atoms.   Reversion  to 
the  ground  state  is  accompanied  by  the  emission  of  visible  light.   Tl 
emission  is  exponential  with  a  decay  constant  of  about  3x]0   seconds. 
The  crystal  is  essentially  transparent  to  its  own  light,  which  is  trans- 
mitted through  the  crystal  to  the  photocathode.   The  photocathode  emits 
electrons  when  illuminated  by  the  light.   These  primary  electrons  are 
accelerated  in  an  electric  field  to  a  metal  dynode  which  produces  second- 
ary  electrons  when  struck  by  the  primary  electrons.   These  secondary 
electrons  are  accelerated  through  another  electric  potential  toward  a 
second  dynode  where  the  process  is  repeated.   There  are  ten  such  dynodes 
in  the  PM  tube,  and  the  result  after  the  last  stage  is  an  electron  mul- 
tiplication of  approximately  10  .    This  electron  current  pulse  is  then 
passed  through  a  load  resistor,  thus  producing  a  voltage  pulse  across  the 
resistor. 

Pulse-height  analysis  is  possible  only  because  the  amount  of  light 
produced  in  the  crystal  is  a  linear  function  of  the  energy  deposited  in 
the  crystal.   The  number  of  primary  electrons  emitted  by  the  photocathode 
is  proportional  to  the  amount  of  light  striking  it.   Consequently,  the 
number  of  electrons  directed  toward  the  first  dynode  is  proportional  to 
the  amount  of  energy  deposited  in  the  crystal.   The  electron  multiplication 
through  the  dynode  stages  is  constant  with  a  stable  high  voltage  power 
supply.   Thus,  the  number  of  electrons  seen  at  the  collector  dynode  is 
still  proportional  to  the  amount  of  energy  deposited  in  the  crystal.   The 
current  pulse  is  allowed  to  pass  through  a  DC  blocking  capacitor  to  a  load 
resistor,  thus  producing  a  voltage  pulse  whose  magnitude  is  proportional 
to  the  energy  deposited  in  the  crystal.   Figure  3  is  a  schematic  diagram  of 
the  PM  tube  and  the  associated  tube  base  electronics  used  in  this  study. 
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FIGURE     3.      PHOTOMULTIPLIER      TUBE     BASE     CIRCUIT. 


The  photon  interactions  and  the  crystal  properties  form  the  basis 
for  interpreting  the  pulse-height  distribution  accumulated  on  the  MCA. 
The  ful]  energy  peak,  or  photopenic,   results  when  all  of  the  energy  of 
the  photon  is  deposited  in  the  crystal  in  a  time  too  short  to  be  resolved 
by  the  detector.   This  may  be  just  one  photoelectric  absorption  or  a 
series  of  Compton  scatters  followed  by  a  photoelectric  absorption. 

The  Compton  distribution  extends  continuously  from  zero  energy  to 
the  Compton  edge,  since  all  energies  between  these  limits  are  possible  in 
the  scattering  process.   The  backscatter  peak  results  from  photons  which 
initially  leave  the  source  in  a  direction  away  from  the  crystal  and  under- 
go a  large  angle  scatter,  about  180  ,  and  are  thus  scattered  back  toward 
the  crystal  and  subsequently  detected.   The  resulting  peak  is  usually  some- 
what broad  because  of  the  cone  of  possible  backscattering  angles  realized 
by  the  source  photons. 

The  characteristic  X-ray  peak  shown  in  Figure  1  is  not  apparent  in 
every  gamma-ray  spectrum,  and  is  not  related  to  the  system  response  to  the 
monoenergetic  photons.   Internal  conversion  processes  which  may  occur  in 
the  isotope  under  study  produce  X  rays  which  may  be  detected  by  the  system. 
These  X  rays  are  produced  by  a  source  which  is  independent  of  the  photon 
source  being  studied.   The  X-ray  peak  must  therefore  be  removed  from  the 
response  curve  associated  with  the  gamma  ray.   The  source  photons  may,  how- 
ever, excite  atoms  in  the  shielding  material  surrounding  the  crystal.   These 
excited  atoms  may  undergo  internal  conversion  processes  and  emit  X  rays 
which  may  be  detected  by  the  crystal.   These  X-ray  peaks  must  remain  in  the 
response  curve  since  they  are  a  consequence  of  the  system  arrangement.   In 
order  to  distinguish  between  these  two  types  of  X  rays   one  must  have  a 
knowledge  of  the  decay  scheme  of  the  isotope  under  study. 
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The  crystal  size  will  also  affect  the  observed  pulse-height  distri- 
bution.  For  instance,  either  one  or  both  of  the  annihilation  photons 
produced  in  pair  production  may  escape  the  crystal  without  further  inter- 
action.  If  either  one  of  the  photons  escape,  a  "one-escape"  peak  will 
appear  in  the  distribution  just  0.511  MeV  below  E  ,  the  source  energy. 
If  both  photons  escape,  a  "two-escape"  peak  will  appear  in  the  distri- 
bution just  1.022  MeV  below  E  .   These  escape  peaks  are  a  consequence  of 
the  fact  that  all  of  the  source  energy,  except  that  which  has  escaped, 
has  been  deposited  in  the  crystal.   As  the  crystal  size  is  increased   the 
probability  for  escape  is  reduced  because  more  material  is  available  for 
the  photon  to  pass  through. 

Another  example  of  the  escape  mechanism  is  the  so  called  iodine 
escape  peak.   When  a  low  energy  gamma  ray  or  X  ray  enters  the  crystal 
it  is  usually  absorbed  very  near  the  surface  of  the  crystal  because  the 
interaction  probability  for  absorption  is  very  high  for  low  energies. 
During  the  absorption  process  iodine  atoms  in  the  detector  are  excited 
and  emit  28  keV  X  rays,  which  have  a  high  probability  for  escape  since 
they  are  produced  very  near  the  detector  surface.   Thus  an  iodine  escape 
peak  appears  in  the  distribution  just  28  keV  below  the  associated  gamma - 
ray  or  X-ray  peak. 

For  a  more  detailed  discussion  of  the  scintillation  process,  the 
reader  is  referred  to  Price  (14) . 

2.2  Photopeak  Parameters 
In  Figure  1  notice  that  the  full  energy  peak,  or  photopeak,  is 
approximately  Gaussian  in  shape.   In  many  cases  the  entire  energy  of  the 
source  photon  is  not  deposited  in  the  crystal  in  one  interaction.   Thus, 
even  though  the  crystal  may  not  be  able  to  resolve  the  various  interactions 
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the  total  amount  of  light  produced  in  the  crysta]  will  vary  about  some 
mean  value.   This  mean  will  correspond  to  the  source  energy,  E  .   Another 
fact  basic  to  the  Gaussian  shape  is  the  statistical  variation  in  the  number 
of  light  photons  produced  by  the  crystal  and  the  number  of  photoelectrons 
produced  by  the  photocathode  for  a  given  amount  of  illumination.  Again,  the 
number  of  photons  and  photoelectrons  produced  will  vary  about  some  mean 
value. 

It  has  been  said  that  the  amount  of  light  produced  in  the  crystal 
is  proportional  to  the  amount  of  energy  deposited.   This  is  not  true  in 
the  strictest  sense,  since  below  about  72  keV  the  crystal  response  is  not 
linear  (13) .   It  may  happen  that  zero  response  occurs  at  some  finite 
energy.   This  non-linearity  causes  the  crystal  to  "add-up"  rapidly 
successive  events  inaccurately;  i.e.,  events  occurring  below  72  keV  do  not 
produce  a  proportional  response.    Consequently,  the  resulting  current 
pulse  produces   a  voltage  pulse  which  may  be  higher  or  lower  than  the  "pro- 
portional" height.   It  has  been  observed  (13)  that  Nal(Tl)  produces  higher 
pulse  heights  than  expected  for  multiple  events.   This  non-linearity  is 
another  reason  why  the  photopeak  is  other  than  a  spike. 

Finally,  the  source  photons  may  undergo  small  angle  scatters  in  air 
or  other  materials  surrounding  the  detector.   These  glancing  collisions 
reduce  the  energy  of  the  photons  somewhat  before  they  enter  the  crystal, 
and  therefore  contribute  to  increased  counts  on  the  lower  energy  side  of 
the  photopeak.   Since  gamma  rays  do  not  gain  energy  in  scattering,  only  the 
lower  energy  side  of  the  peak  is  affected  and  the  distribution  is  not  exact- 
ly symmetric,  although  in  most  cases  the   asymmetry  is  not  noticeable. 

As  previously  noted,  the  data  are  accumulated  on  a  multichannel 
differential  pulse-height  analyzer  (MCA) .   The  baseline  control  is  set  such 
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that  zero  energy  corresponds  approximately  to  channel  one.   The  gain  con- 
trol is  set  such  that  the  energy  range  of  the  MCA  extends  beyond  the 
source  photon  energy.   The  accumulated  energy  spectra  are  then  presented 

in  discrete  form  as  number  of  counts,  Y.,  in  channel  i,  where 

1 

i  =  1,  2 ,...,400.   The  data  are  thus  represented  as  a  histogram.   The 
histogram  displays  the  count  rate  as  a  function  of  pulse  height  (energy) . 
It  is  important  to  remember  that  each  channel  in  the  histogram  has  a 
finite  energy  width,  AE.   If  channel  i  is  representative  of  energy  E.,  then 
all  detected  photons  with  energies  between  E.  +  AE/2  are  registered  in 
channel  i.   Notice  that  the  center  of  the  channel  is  chosen  as  represent- 
ative of  all  photon  energies  registered  in  that  channel.   In  other  words, 
all  photons  with  energies  between  E.  +  AE/2,  which  are  detected,  are 
assumed  to  have  an  energy  E.. 

Since  the  photon  interaction  cross  sections  in  the  crystal  are  con- 
tinuous functions  of  energy,  the  pulse-height  distribution  is  a  continuous 
function.   If  Y(x)  represents  the  pulse-height  distribution,  where  the  x 
are  the  pulse-height  units,  the  count  rate,  Y.,  in  channel  i  is  given  as 


,    Ax 
x.   +  -x— 

f     1 

Y.    = 

X 

Y(x)dx 

, 

Ax 
xi        2 

(2.2) 


where  x.  is  the  mean  pulse  height  in  channel  i.   It  is  assumed  that  Ax  is 
sufficiently  small  such  that 


Y.  =  Y(x.)Ax.  (2.3) 

l      l 

As  previously  discussed,  the  photopeak  may  be  represented  as  a 
Gaussian  function, 


13 


t   -      )2 
Y(x)  -  Y(x  )  exp[  -^--SlL     j  m  (2./,) 

°         Bo 

Y(x  )  is  the  maximum  count  rate  in  the  photo-peak;  x   is  the  pulse  height 
o  o 

at  which  the  maximum  occurs;  this  is  not  to  be  confused  with  a  channel 

number.   W  is  the  full  width  of  the  photopenic  at  half  maximum  (i'l.'HM)  as 
o 

given  by 

W  -  2(B  ln2)1/2   .  (2.5) 

o      ° 

Equation  (2.4)  may  be  used  to  calculate  individual  points  on  a  Gaussian 
distribution  by  rewriting  it  as 


Y(x.)  =  Y(x  )  exp  [  -^L_x-^  ]   .  (2.6) 

l        o  BQ 


Using  the  approximation  shown  in  equation  (2.3)  in  conjunction  with 
equation  (2.6)  the  photopeak  data  may  be  represented  as 

-(x.-x  y 

Yi  *  Yo  £XP  [    B      ]   '  (2*7) 

o 

where  Y  =  AxY(x  ) . 
o        o 

Equation  (2.7)  is  suitable  for  use  in  a  weighted  polynomial  least 
squares  analysis  if  one  linearizes  the  equation  by  taking  the  natural 
logarithm  of  both  sides.   The  result  is  an  equation  of  the  form 

ax2  +  bx  +  c  =  f(x)  (2.8) 
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where, 

-1 

a  =  B~ 
o 

b  --^o 

B 

o 

x  2 
c  =  lnY  -  Xo    .  (2.9) 

°   B~ 
o 

Solving  equation  (2.8)  for  a,  b,  and  c  will  yield  the  desired  fitting 
parameters  x  ,  Y  ,  and  B   after  some  simple  manipulations. 

The  data  in  each  channel  will  vary  statistically,  and  it  is 
assumed  that  the  statistical  distribution  function  of  each  channel  is 
Gaussian.   The  weighting  factors  used  in  the  least  squares  analysis  are 
then  given  as  the  inverse  of  the  variance  of  f(x),  where  f(x)  =  lnY(x). 
Standard  propagation  of  error  techniques,  as  discussed  by  Beers  (2), 
yield  the  weighting  factor  in  the  log  domain  as: 

W(x)   =  Y(x).  (2.10) 

The  pulse-height  resolution  of  the  spectrometer  system  is  an  over- 
all measure  of  the  system  performance  (3,  16).   The  various  factors  which 
contribute  to  the  observed  resolution  are  discussed  by  many  authors, 
(14,3,17,8,4,16)  and  are  only  briefly  stated  here.   The  factors  may  be 
categorized  in  four  main  areas,  (i)  variance  in  the  emission  of  light 
photons  by  the  scintillator,  (ii)  variance  in  the   collection  of  the  light 
photons  by  the  photocathode,  (iii)  variance  in  the  emission  of  electrons 
from  the  photocathode  and  collection  at  the  first  dynode,  and  (iv)  variance 
in  the  electron  multiplication  process  through  the  dynode  stages.   All  of 
these  factors  are  statistical  in  nature  and  may  be  represented  mathemati- 
cally (3,4,17).   All  of  these  factors  operate  in  a  staged  process  to  dis- 
tort the  photopeak  response,  the  result  being  a  Gaussian  shape. 
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B  ,  the  width  parameter,  is  used  to  calculate  the-  resolution  of 
o 

the  system.   Equation  (2.5)  allows  the  calculation  of  W  ,  which  may  be 
stated  in  units  of  energy  or  channels.  X    ,  the  photopeak  pulse  height, 
also  corresponds  to  an  energy,  i.e.,  E  ,  the  source  energy.   The  system 
resolution,  R,  is  then  stated  as  a  percent  which  is  calculated  as  shown 
in  equation  (2.11). 

R  =  ^J°   .  (2.11) 

x 

o 

A  perfect  detector  would  display  only  a  spike,  with  a  resolution  of  zero. 
In  most  cases,  a  resolution  of  about  7  percent  for  Cs-137  is  considered 
good  for  Nal(Tl). 

2.3   Spectrum  Analysis 

The  first  step  in  the  data  analysis  is  the  least  squares  fit  of  all 
the  raw  data  to  determine  the  peak  pulse  height,  x  ,  as  well  as  the  other 
basic  parameters  previously  mentioned.   Also,  composite  Co-60-Cs-Ba-137 
quality  control  spectra  must  be  fitted  in  order  to  determine  the  zero  en- 
ergy channel  of  the  calibration  spectra.   This  is  a  very  important  step, 
because  without  knowledge  of  where  the  zero  intercept  is  for  the  calibra- 
tion spectra,  the  data  accumulated  may  well  be  useless.   The  quality  con- 
trol spectra  are  fitted  in  the  same  manner  as  the  experimental  spectra.  In 
addition,  a  least  squares  straight  line  is  computed  to  relate  x  and  E  , 
thus  determining  the  calibration  slope,  keV/channel,  and  hence  the  zero-en- 
ergy intercept. 

Once  x  and  the  zero  energy  intercept  are  determined  for  the  cali- 
bration spectra,  it  is  necessary  to  shift  the  data  in  a  linear  manner  such 
that  all  of  the  data  are  adjusted  to  the  same  gain  and  base  line  settings. 
A  very  small  gain  shift  may  not  change  the  spectral  data  appreciably,  but 
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the  effect  on  the  spectrum  statistics  is  quite  noticeable.   Consequently 
a  gain  shift  of  even  0.5  channels  must  be  accounted  for  (5). 

In  this  study  nine  different  isotopes  were  used  to  measure  re- 
sponse functions  at  energies  between  46  keV  and  1.114  MeV.   Several  sets 
of  data  were  taken  for  each  collimator  aperture  and  for  each  position. 
The  experimental  arrangement  will  be  discussed  in  more  detail  in  the 
following  sections.   A  total  of  approximately  750  spectra  were  obtained 
over  an  eight  month  period.   Consequently,  all  of  the  spectra  required 
shifting,  since  gain  drifts  and  shifts  do  occur  daily. 

As  previously  stated,  the  shifts  in  the  measured  spectra  are  linear 

because  it  is  assumed  that  the  MCA  operates  in  a  linear  manner.   If  x..  is 

the  photopeak  pulse  height  in  the  measured  spectrum  and  x?  is  the 

corresponding  pulse  height  in  the  adjusted  spectrum,  the  gain  shift  ratio, 

r,  may  be  given  by 

x 
r  =  —  (2.12) 

Xl 


assuming  that  there  is  no  zero  energy  shift.  Since  the  variables  x  and 
x~  are  continuous  functions  of  energy,  it  is  possible  to  calculate  the 
new  count,  y(x„),  at  x_  as 

yOO 

y(x2)  =  — ^-  .  (2.13) 

The  gain  shift  ratio  used  maintains  the  total  area  under  the.  curve;  i.e., 
the  total  number  of  counts  remains  the  same  in  the  two  spectra  (8) . 

For  the  purposes  of  illustration,  assume  that  two  measured  pulse-height 
distributions  Y(x)  and  Y'(x),  are  available  for  the  same  experiment.   How- 
ever, let  the  two  distributions  be  on  different  pulse-height  scales.   Let 
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Y(x)  be  represented  by  m  channels  of  data  and  V'(x)  by  n  channels  of  data, 

with  m  less  than  n.   Both  distributions  cover  the  same  energy  range,  but 

the  channel  width  for  Y'(x)  is  smaller  than  that  for  Y(x).    In  order   to 

compare  the  two  distributions  directly,  the.  channel  widths  of   each  must 

be  the  same.   Since  the  energy  range  is  the  same,  J-  channels  in  the  Y'(x) 

distribution  must  be  equivalent  to  one  channel  in  the  Y(x)  distribution. 

The  shifting  ratio,  r,  is  — .   In  order  to  make  the  pulse-height 

m 

scales  for  the  two  distributions  the  same,  Y'(x)  is  divided  by  r  and  the 
channel  numbers  for  Y'(x)  are  multiplied  by  r.   Even   though  the  pulse- 
height  scales  are  the  same  after  the  shift,  the  data  in  the  shifted  spec- 
trum do  not  occur  at  the  same  pulse  heights  as  the  data  for  Y(x).   Conse- 
quently, an  interpolation  procedure  must  be  introduced  in  order  to  calcu- 
late the  count  rate  at  the  proper  pulse  heights.   The  interpolation  is 
accomplished  by  choosing  a  data  point  nearest  the  proper  pulse-height  and 
a  point  on  either  side.   These  three  points  are  used  to  calculate  the  con- 
stants in  a  quadratic  equation.   This  equation  is  then  used  to  calculate 
the  count  rate  at  the  required  pulse  height.   The  process  is  repeated  m 
times  until  both  distributions  match  over  the  entire  pulse-height  range. 
To  change  the  zero  energy  intercept   of  the  measured   spectrum  as 
well  as  the  gain,  equation  (2.12)  becomes 

r-oJir  (2-14) 

Where  6  is  the  desired  shift  in  the  zero  energy  intercept  of  the  measured 
spectrum.   Equation  (2.13)  remains  unchanged. 

2.4   Simulation  of  Infinite  Plane  Sources  Using   Point  Sources 

Point  source  spectra  were  accumulated  for  sources  placed  at  equally 


18 


spaced  intervals  off  the  center  line  of  the  collimated  detector.   A 
schematic  diagram  of  the  geometry  is  shown  in  Figure  4.   For  the  two-inch 
diameter  collimator   the  sources  were  placed  at  1/2"  intervals  out  to  a 
radius  of  4  1/2"  to  5".   The  intervals  were  1/4"  out  to  a  radius  of  1"  to 
1  1/4"  for  the  one-half  inch  diameter  collimator.   The  last  radius  point 
used  was  the  position  at  which  the  count  rate  in  the  photopeak  was  only  2 
or  3  counts  per  minute  above  background.   Thus,  effectively,  the  contri- 
bution to  the  response  function  out  to  an  infinite  distance  was  measured. 
Figure  5  is  a  typical  representation  of  the  photopeak  amplitude,  Y  ,  as 
a  function  of  the  distance  off  the  detector  center  line.   Notice  that  the 
amplitude  in  the  last  position  is  very  close  to  zero  and  that  it  is  at 
least  two  orders  of  magnitude  below  the  count  rate  measured  at  the  center 
position. 

In  order  to  represent  the  infinite  plane  response,  the  point  source 
data  described  above  must  be  numerically  integrated  in  order  ,  by  superpo- 
sition,  to  generate  the  response  function  for  a  disk  source  of  effectively 
infinite  radius.   Before  the  data  can  be  integrated,  however,  they  must  be 
corrected  for  radioactive  decay  and  background,  and  shifted  so  that  the 
photopeaks  for  one  particular  isotope  all  appear  at  the  same  pulse  height, 
and  also  shifted  so  that  all  of  the  spectra  have  the  same  zero  energy  in- 
tercept.  Finally,  each  spectrum  must  be  normalized  to  the  experimental 
isotope   source  strength.   All  of  these  corrections  and  normalizations  are 
accomplished  by  a  code  written  for  that  purpose  called  SPECDATA. 

The  numerical  integration  of  the  corrected  normalized  spectra  is  per- 
formed by  DISK  which  utilizes  the  WATE  subroutine  (10).   WATE  uses  Simpson's 
rule  for  an  odd  number  of  integration  points  and  a  modified  Simpson's  rule 
for  an  even  number  of  points  to  construct  a  list  of  integration  weights  for 


19 


a, 

•H 

JZ 
03 

c 
o 

•H 
•U 
nJ 

rH 

0) 

c 

CO 

rH 

a- 

cu 
o 
u 

o 

CO 

I 

u 
o 


o 
o 

4-1 

o 

(3 

o 

•H 


<u 

CO 
CU 
V-i 

ex 

CU 

u 


u 

CO 

§ 

a 

CO 


cu 

a 

bO 
•H 


20 


cu 

Pi 


C 

o 
u 

0) 

> 

•H 

-u 
cd 

rH 
CU 

Pi 


Figure   5. 
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Distance  off  center  line,  (inches) 

Typical  photopeak  response  as  a  function  of  the  distance  off 
center,  Mn-54,  half-inch  collimator. 
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equally  spaced  data.   When  only  two  points  are  available  for  the  in- 
tegration, a  trapezoidal  approximation  is  applied.   Equation  (2.15)  shows 
the  mathematical  scheme  used  for  the  numerical  integration. 
'2tt    rR 
d(Ji 


N 


:.(r)rdr  =  2tt   Y   f.  (r.)r.W.  ;  j  =  1,2,. ..,400.    (2.15) 
1  ,L.     j   i  i  ±   J         ' 

J  i=l 


N  is  the  number  of  radial  points  used  in  the  integration,  providing  N-l 
intervals.   W.  are  the  Simpson  integration  weights.   The  integration 
weights  are  returned  from  WATE  to  MAIN  in  DISK  and  the  data  in  each  chan- 
nel of  each  spectrum  are  multiplied  by  the  appropriate  weight.   The  multi- 
plied spectra  are  then  added  together  channel  by  channel  to  yield  the 
system  response  to  a  radial  line  source.   The  system  response  is  azimuthally 
symmetric,  consequently  multiplying  the  line  source  response  by  2tt  yields 
the  system  response  to  a  disk  source  of  effectively  infinite  radius. 

After  the  integration  of  the  data  for  a  particular  isotope  is  com- 
pleted, the  photopeak  of  the  response  function  is  fitted  to  a  Gaussian, 
for  a  final  check,  in  order  to  make  sure  that  the  integration  process  did 
not  change  the  position  of  the  photopeak  appreciably. 


22 


3.0   EXPERIMENTAL  PROCEDURE 

3.1  Facilities 

The  goal  of  the  experiment  was  to  determine  the  response  of  the 
system  to  gamma  rays  from  an  infinite  plane  monoenergetic  source.   As  pre- 
viously discussed,  the  infinite  plane  was  simulated  by  using  point  sources 
in  a  plane  about  100  centimeters  from  the  crystal  face.   The  crystal  was 
placed  inside  a  cylindrical  lead  collimator.   Figure  6  is  a  machine  draw- 
ing of  the  collimator  assembly  with  the  appropriate  details.   Such  a  de- 
vice is  necessary  for  many  measurements  where  angular  and  energy  distri- 
butions are  of  interest.   The  collimator  has  two  degrees  of  freedom;  the 
base  will  rotate  through  360  ,  while  the  barrel  itself  will  rotate  through 
90  .   There  are  three  collimator  plugs  available  with  the  assembly,  a  solid 
lead  blank  for  use  in  background  measurements,  a  one-half  inch  diameter 
collimator,  and  a  two-inch  diameter  collimator.   The  half-inch  collimator 
has  an  angular  resolution  of  1.85   as  measured  from  the  center  of  the  de- 
tector face.   The  two-inch  collimator  has  a  resolution  of  7.39  . 

Figure  7  is  a  block  diagram  of  the  experimental  arrangement  and  the 
equipment   used.   The  infinite  plane  was  simulated  by  placing  the  point 
sources  on  a  wooden  rack  which  was  parallel  to  the  collimator  face.   Fig- 
ure 8  is  a  photograph  of  the  collimator  assembly  and  source  rack.   The 
meter  stick  seen  on  the  left  was  used  to  measure  the  source  position  rel- 
ative to  the  center   line  of  the  collimator   entrance.   A  plumb  bob  was 
hung  from  the  top  of  the  rack  such  that  the  plumb  bob  string  intersected 
the  horizontal  string  at  the  center  of  the  collimator  entrance  as  project- 
ed on  the  source  plane  described  by  the  source  rack.   The  center  position 
was  checked  and  calibrated  at  least  twice  daily  with  an  optical  level  to 
insure  reproducibility  of  the  results. 
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1.  Pb  collimator 

2.  PM  tube 

3.  3"x3"  Nal(Tl)  crystal 

4.  HV  power  supply 

5.  Pre-amplif ier 

6.  Pre-amplif ier  power  supply 

7.  Input-output  control  unit 

8.  400  channel  pulse-height 
analyzer 

9.  Point  source  locations 
10.  Optical  level 


Figure  7.   Block  diagram  of  equipment  set-up. 
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Figure  8.   Collimator  assembly  and  source  rack. 
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The  potentiometer  shown  in  the  lower  left  of  Figure  8  was  used  in 
conjunction  with  a  copper-constantan  thermocouple  to  monitor  the  crystal 
environment  temperature  inside  the  collimator.   It  is  desirable  to  main- 
tain a  constant  crystal  and  PM  tube  temperature  in  order  to  minimize  gain 
drifts  and  shifts.   Also,  the  thermionic  emission  current  from  the  dynodes 
and  photocathode  should  be  kept  as  low  as  possible.   Price  (14)  indicates 
that  the  thermionic  current  may  be  reduced  by  about  a  factor  of  two  for 
every  ten  degree  centigrade  drop  in  tube  temperature.   Since  it  is  not 
convenient  to  cool  the  collimator,  it  was  decided  to  maintain  a  tempera- 
ture as  near  room  temperature  as  possible  inside  the  collimator.   This 
was  accomplished  with  a  constant  temperature  water  bath  heater.   Heated 
water  was  pumped  directly  into  the  collimator  through  copper  tubes  im- 
bedded in  the  lead  (see  Figure  6) . 

The  detector  used  was  a  Harshaw  Integral  Line  Assembly,  type  12S12, 
consisting  of  a  3"  x  3"  Nal(Tl)  crystal,  PM  tube  and  associated  tube  base 
electronics.   The  PM  tube  was  powered  by  a  John  Fluke,  model  412A,  high 
voltage  power  supply,  operated  at  +1100  Volts.   The  signal  from  the  PM 
tube  was  passed  through  a  Hamner  model  N-356  cathode  follower  pre-amplif ier , 
powered  by  a  Kepco  model  HB2AM  voltage  regulated  power  supply  operated  at 
+250  Volts.   The  output  of  the  pre-amplif ier  was  fed  through  approximately 
eighty  feet  of  RG-62U  signal  cable  to  the  input  terminal  of  a  TMC  model 
402-6  400  channel  MCA.   Output  from  the  MCA  was  controlled  by  a  TMC  model 
520  type-punch-read  control  unit,  which  permitted  read-out  by  an  IBM  type- 
writer or  a  Tally  paper  tape  perforator,  and  read-in  by  a  Tally  paper  tape 
reader. 

3.2   Source  Preparation 

In  this  study  nine  isotopes  were  used  as  sources  for  monoenergetic 
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ganuna  rays;  lead-210  (0.046  MeV) ,  cadraium-109  (0.088  MeV),  cerium-139 
(0.166  MeV),  mercury-203  (0.279  MeV),  indium-tin-113  (0.392  MeV), 
strontium-85  (0.513  MeV),  cesium-barium-137  (0.662  MeV) ,  manganese-54 

(0.840  MeV),  and  zinc-65  (1.114  MeV).   All  of  the  isotopes  were  purchased 
commercially  in  liquid  form  and  were  certified  to  be  99+  percent   pure. 
Lead-210  was  purchased  from  New  England  Nuclear   Corp.,  575  Albany  Street, 
Boston,  Mass.;  all  other  isotopes  were  purchased  from  Nuclear   Science 
and  Engineering  Corp.,  P.  0.  Box  10901,  Pittsburgh,  Pa. 

Before  the  sources  could  be  used  in  the  experiments   it  was 
necessary  to  transform  them  into  "point"  form.   Source  mounts  were  con- 
structed from  1/16"  sheets  of  ordinary  plastic  and  consisted  of  a  1"  O.D. , 

3/4"  I.D.  circular  section  mounted  on  a  two  inch  long  handle  to  facilitate 

2 
handling.   A  thin  sheet  of  transparent  mylar  film, 3. 5  mg/cm  ,  was  sealed 

to  one  side  of  the  planchet  with  Walther's  "Goo"  (manufactured  by  Walthers 

Specialty  Inc.,   Milwaukee,  Wisconsin)  which  was  found  to  give  a  strong, 

leak-proof  seal.   Typically,  the  planchet  weight  before  loading  was  about 

0.48  grams.   The  procedure  was  then  to  draw  the  source  solution,  usually 

less  than  0.1  ml,  from  the  source  bottle  by  means  of  a  micrometer  buret 

(manufactured  by  RGI  Inc.,  Box  732,  Vineland,  New  Jersey,  model  S1200A) , 

with  a  0.2  ml  capacity  and  0.0002  ml  divisions.   Figure  9  shows  a  typical 

source  bottle,  unloaded  planchet,  and  the  disassembled  buret. 

The  entire  solution  was  drawn  from  the  source  bottle  and  the  volume 

recorded.   The  buret  tip  was  then  placed  above  the  center  of  the  planchet 

and  clamped  into  position.   Approximately  1/20  of  the  total  volume  was 

placed  on  the  planchet  one  drop  at  a  time,  the  buret  was  then  rotated  away 

from  the  planchet  and  an  infra-red  lamp  placed  in  position  over  the  source 

solution,  which  was  allowed  to  evaporate,  under  a  hood,  until  just  the  source 
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residue  was  left  on  the  mylar  film.   This  process  was  repeated  until  the 
desired  amount  of  material  was  dried  onto  the  planchet. 

The.  planchet  was  then  removed  to  a  counting  cavity,  which  is  de- 
scribed later,  and  counted.   After  counting,  the  entire  process  was  repeat- 
ed and  the  source  was  counted  again.   Then  the  remaining  solution  was  dried 
onto  the  planchet,  one  drop  at  a  time.   It  is  important  that  this  process 
be  completed  in  one  operation,  so  that  the  solution  in  the  buret  does  not 
have  an  opportunity  to  evaporate.   An  H  SO,  +  KMnO,  solution  was  used  to 
clean  the  glass  part  of  the  buret  immediately  after  the  completion  of  a 
source  preparation.   This  insured  that  the  buret  did  not  become  radio- 
actively  contaminated  in  any  measurable  manner. 

3.3  Procedure 

The  analyzer  equipment  and  high  voltage  power  supply  were  operated 
by  110  Volt,  60   Hz  power.   An  intermatic  switch  automatically  turned  the 
equipment  on  at  least  two  hours  before  the  experiment  day  started.    This 
allowed  the  equipment  to  warm  up  and  the  high  voltage  power  supply  to 
stabilize.   Another  intermatic  switch  turned  on  an  air  conditioner  at  about 
the  same  time  so  the  room  remained  cool  and  the  equipment  did  not  become 
over  heated.   The  pre-amplif ier  power  supply  was  left  on  24  hours  a  day 
since  no  intermatic  switching  was  available. 

At  the  beginning  of  the  experiment  day  a  series  of  standard  tests 
were  performed  on  the  MCA  (19)  to  determine  if  it  was  operating  properly. 
In  addition,  a  quality  control  run  using  a  cobalt-60  and  a  cesium-barium-137 
source  was  performed  in  order  to  determine  the  zero  energy  intercept.   Such 
a  run  was  performed  about  every  two  hours  in  order  to  maintain  an  accurate 
estimate  of  the  zero  energy  intercept.   The  gain  of  the  MCA  was  set  so  that 
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the  photopeak  of  the  calibration  source  under  study  was  in  the  highest 
possible  channel  with  the  restrictions  that  the  noise  in  the  first  few 
channels  were  maintained  at  a  minimum  and  that  no  high  energy  counts 
were  lost  beyond  channel  400. 

The  next  step  in  the  procedure  was  the  positioning  of  the  colli- 
mator and  the  source  rack.   This  was  one  of  the  more  critical  steps  in 
the  experiment,  since  a  relatively  small  error  in  the  positioning  could 
produce  a  large  error  in  the  accumulated  spectra  and  hence  in  the  final 
response  function.   The  collimator  was  checked  each  morning  to  be  sure 
that  the  barrel  was  level  and  that  the  collimator  face  was  parallel  to 
the  source  rack.   The  center  position  was  checked  and  located  with  the 
aid  of  an  optical  level  mounted  on  a  tripod  located  about  15  feet  behind 
the  source  rack.   The  source  rack  was  also  checked  to  be  sure  that  it  was 
level.   This  entire  process  was  repeated  at  least  twice  daily.   The 
positioning  procedure  allowed  source  placement  to  within  1/8   inch,  the 
diameter  of  the  largest  source  dot. 

At  least  two  spectra  were  accumulated  at  each  source  position  for 
each  calibration  source  used.   A  minimum  live  time  was  chosen  so  that  at 
least  1,000  counts  were  registered  in  the  maximum  channel  of  the  photopeak. 
This  required  runs  of  about  10  minutes  at  the  center  position,  and  runs  as 
long  as  1,000  minutes  at  the  farthest  radial  position.   Background  runs 
were  usually  taken  during  the  night  for  about  1,000  minutes.   These  runs 
were  preceded  and  followed  by  quality  control  runs. 

Five  of  the  nine  sources  used,  Pb-210,  Hg-203,  In-Sn-113,  Cs-Ba-137, 
and  Zn-65,  emitted  electrons.   These  electrons  could  be  detected  and  con- 
tribute to  the  response  which,  of  course,  is  undesirable.    In  order   to 
eliminate,  these  effects  from  the  response,  a  1/4"  thick  lucite  beta  shield 
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was  placed  in  front  of  the  crystal  face  to  absorb  these  electrons.   How- 
ever, the  shield  does  have  an  effect  on  the  gamma-ray  flux  incident  on 
the  shield.   The  total  spectrum  must  be  corrected  for  the  flux  attenuation, 
These  corrections  are  discussed  in  detail  in  Appendix  B. 
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4.0  RESULTS 


4.1  Pvesponse  Function  Normalization 
The  response  functions  must  be  normalized  in  some  consistent  manner 
so  they  will  be  useful  to  other  experimenters.   In  this  study  the  responses 
are  normalized  to  the  unit  uncollided  flux  per  steradian  at  the  source 

!in  a  direction  perpendicular  to  the  source  plane.   The  source  strength  of 

2     -1 
the  plane  is  given  as  S.  (cm  -sec)   .   The  uncollided  flux  is  then  given 

as 

g 
,uc    A  ,   2       .-1  .   _N 

<fc   =  7—  (cm  sec  sr)     .  (4.1) 

4tt 

Before  the  point  source  responses  were  numerically  integrated  to 
give  the  infinite  plane  response  functions,  they  were  normalized  to  a  unit 
point  source  strength,  which  had  been  previously  determined  by  experiment. 
The  integration  process  itself  inherently  normalizes  the  plane  source  re- 
sponses to  a  unit  area  source  strength.   Therefore,  in  order  to  complete 

the  normalization,  the  response  functions  are  merely  multiplied  by  4tt  .   The 

2        -1 
units  of  the  resulting  response  functions  are  counts/sec  per  (cm  sec  sr) 

All  of  the  response  functions  presented  are  on  an  energy  scale  of  3.5  keV  per 

channel. 

4.2   Two-Inch  Collimator  Results 
Figure  10  shows  the  response  function  obtained  for  1.114  MeV  gamma  rays 
This  result  is  considered  typical  of  all  nine  two-inch  response  functions  in 
that  the  data  are  smooth,  which  indicates  good  statistical  results.   Heath  (7) 
shows  a  peak-to-valley  ratio  of  approximately  28  for  a  point  source  response 
function  of  the  same  energy  obtained  in  non-collimated  geometry.   Figure  10 
shows  a  similar  ratio  of  approximately  16.   The  difference  is  due  to  colli- 
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mator  effects  such  as  small  angle  scattering  in  the  collimator  barrel  and 
the  lucite  beta  shield.   These  effects  tend   to  fill  the  valley  between 
the  photopeak  and  the  Compton  edge.   The  only  error,  however,  is  caused 
by  scattering  which  occurs  in  the  beta  shield.   The  other  effects  intro- 
duced by  the  collimator  are  part  of  the  system  response.   This  error,  how- 
ever, is  no  larger  than  1  percent  in  Figure  10. 

The  spectrum  shape  and  other  characteristics  are  in  good  agreement 
with  those  observed  by  Heath  in  his  measurements,  except  for  the  following 
two  facts.   Zn-65  was  used  as  the  source,  and  this  particular  isotope 
emits  a  positron  which  is  detected  by  the  crystal.   When  the  positron  re- 
combines  with  an  electron  in  the  crystal  two  0.511  MeV  photons  are  created. 
Since  the  positron  is  not  very  energetic,  recombination  takes  place  near 
the  crystal  face  and  one  of  the  photons  thereby  created  has  a  very  good 
chance  of  escaping  the  detector.   The  other  photon,  however,  has  a  very 
good  chance  of  being  detected;  consequently  a  0.511  MeV  annhilation  peak 
appears  in  the  spectrum.   This  peak  is  not  part  of  the  system  response  and 
must  therefore  be  removed  from  the  spectrum.   After  the  removal  the  re- 
sponse is  no  longer  the  response  to  a  Zn-65  source,  but  rather  the  system 
response  to  a  source  of  1.114  MeV  photons.    The  second  characteristic  to 
be  noted  is  the  small  peak  near  channel  20.   This  peak  occurs  at  approxi- 
mately 70  keV  on  the  energy  scale  and  is  caused  by  lead  X  rays.   The 
1.114  MeV  photons  interact  with  the  collimator,  which  is  constructed  pri- 
marily of  lead.   The  source  photons  are  energetic  enough  to  eject  a  K-shell 
electron  from  a  lead  atom  and  the  subsequent  internal  conversion  process  of 
the  atom  causes  the  emission  of  a  72  keV  X  ray  which  is  detected  by  the 
crystal.   This  peak  is  apparent  in  nearly  all  of  the  response  functions 
and  is  a  part  of  the  system  response.   The  peak  does  not  appear  at  exactly 
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72  keV  on  the  energy  scale  because  of  crystal  non-linearities  at  and  below 
this  energy. 

Figure  11  is  a  composite  drawing  of  all  nine  two-inch  collimator 
response  functions.   The  general  trends  shown  are  relatively  obvious.   The 
photopeak  amplitudes  decrease  with  increasing  energy,  which  is  expected 
since  the  absorption  cross  section  decreases  with  energy.   This  is  shown 
more  definitely  in  Figure  12.   The  lead  X-ray  peaks  are  apparent  in  most 
of  the  spectra  around  channel  15.   These  peaks  are  not  obvious  in  the  two 
lower  energy  response  functions,  either  because  they  are  not  produced  by 
the  low  energy  source  photons  or  because  they  are  masked  by  poorly  resolved 
iodine  escape  peaks. 

Figure  12  shows  the  photopeak  heights  for  the  two-inch  collimator  as 
a  function  of  energy.   At  lower  energies  photo   absorption  is  the  dominant 
photon  interaction  process  in  the  crystal  detector  and  increases  very 
rapidly  with  decreasing  energy.   In  the  middle  energy  range  Compton 
scattering  becomes  the  dominant  process  and  the  absorption  process  decreas- 
es exponentially  in  Nal(Tl).   Figure  12  is  also  a  good  indication  of  the 
consistency  of  the  two-inch  collimator  data. 

The  width  parameter  of  the  photopeaks,  B  ,  is  shown  in  Figure  13  as 
a  function  of  energy.   Above  approximately  300  keV,  B  seems  to  be  increas- 
ing in  a  linear  manner,  indicating  that  the  system  resolution  is  improving 
with  energy.   As  energy  decreases,  the  resolution  decreases;  i.e.,  the 
system  cannot  resolve  peaks  which  are  closely  spaced.   Inspection  of  Figure 
13  indicates  that  B   is  tending  toward  zero  at  zero  energy,  which  means 
that  the  W  is  zero  at  zero  energy. 

The  actual  response  functions  are  tabulated  in  Appendix  D  with  their 
appropriate  variances.   The  standard  deviations  of  the  response  functions 


36 


_|    UJ 

°! 

i 

O   UJ 

z  i- 


U 

z 
z 
< 

X 

o 

to  ^ 

z  ■>-- 

o  m 

H  *> 

o  « 

-  UJ  CD 

P  _J  o 

U.  <    (y> 

o  — 

UJ  CO 

co       cr 

z  >■  y 

0.EO 

co  uj  t- 

1x3  Z   <-> 

t    l'l     O 


3ivu  iNnoo  3Aiivnay 


37 


10.0 


9.0 


8.0 


7.0 


tJ   6.0 

Pi 


5.0 


4.0 


3.0 


2.0 


1.0 


.J 


0.2 


0.4 


0.6      0.8       1.0 
Energy,  (MeV) 


1.2 


Figure  12.   Photopeak  amplitude  as  a  function  of  energy,  two-inch  collimator, 
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Figure  13.   B   as  a  function  of  energy  for  two-inch  collimator 
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for  the  two-inch  collimator  average  somewhat  below  5  percent.   These 
standard  deviations  include  the  error  in  the  source  calibrations  and  the 
random  error  associated  with  the  counting  process.   The  standard  deviations 
are  indicative  of  the  statistical  processes  of  detection  only  and  no  est  i- 
mate  of  the  systematic  error  is  included.   These  systematic  errors  could 
be  caused  by  errors  in  the  source  positioning  on  the  source  rack,  com- 
puter round-off  error  in  the  data  analysis,  and  errors  that  may  have 
occurred  in  measuring  the  source  to  detector  distances.   It  is  believed, 
however,  that  these  types  of  errors  are  small  in  the  two-inch  collimator 
experiments  and  that  the  standard  deviation  would  not  be  increased  by  more 
than  one  or  two  percent  if  these  errors  could  be  included. 

4.3   One-half  Inch  Collimator  Results 

The  half-inch  response  functions  were  accumulated  in  the  same  manner 
as  the  two-inch  responses  and  normalized  in  the  same  manner  to  the  unit 
uncollided  flux  per  steradian  and  are  presented  on  the  same  energy  scale 
as  the  two-inch  responses,  i.e.,  3.5  keV/channel.   There  are  only  eight 
response  functions  however,  instead  of  the  nine  obtained  for  the  two-inch 
collimator  experiments.  The  Hg-203  response  function  was  not  obtainable 
for  the  half-inch  experiments  because  of  MCA  malfunctions  which  invali- 
dated most  of  the  data  obtained. 

If  a  point  source  is  placed  on  the  collimator  center  line,  the  ob- 
served countrate   for  the  half-inch  collimator  should  be  approximately  a 
factor  of  16  lower  than  that  observed  with  the  two-inch  collimator  because 
the  solid  angle  subtended  is  factor  of  16  smaller. 

It   can  be  shown  however,  that  the  half-inch  response  to  a  plane 
source  is  approximately  a  factor  of  256  less  than  the  two-inch  collimator 
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response  to  a  plane  source.   Let  R..  and  R  represent  the  two-inch  and  half- 
inch  responses,  respectively,  for  a  plane  source,  and  C  and  C   represent 

the  responses  for  a  point  source. 

'r. 


R±   =  2tt 


C  (r)rdr 


R2  =  2tt 


(4.2) 


C2(r)rdr 


where  r..  and  r„  are  the  radii  of  the  plane  sources  used  for  the  two-inch 
and  the  half-inch  collimators.   For  the  purposes  of  illustration,  it  is 
assumed  that  the  C.  are  independent  of  r.    In  these  experiments  r..  was 
approximately  four  times  r_.   We  then  have  the  following  relationships, 


Rl  *  "1  Cl 


R2  =  ^r2  C2 


rl  =  4r2 


C1  =  16 C2 


(4.3) 


The  ratios  of  the  responses  are  computed  as  shown  below; 
R, 


r  2C 
rl  Cl 


"^2 


(16r2  ><16C2) 
^7 


=  256 


(4.4) 


It  was  experimentally  observed  that  the  corresponding  point  source  re- 
sponses for  the  two  collimators  were  related  by  approximately  a  factor  of 
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sixteen. 

Half-inch  collimator  experiments  should  also  show  improved  system 
resolution  because  scattering  in  the  collimator  barrel  should  be  reduced 
because  of  less  surface  area  available  for  the  scattering.   Improved 
resolution  was  observed  in  this  study  as  is  shown  in  Figure  14.   The  half- 
inch  resolution  seems  to  be  improved  by  approximately  13  percent  over  the 
two-inch  experiments.   It  should  be  noted  that  the  percent  resolution  of 
both  collimators  are  smooth  functions  in  the  energy  range  of  interest. 

Figure  15  represents  the  photopeak  amplitude  as  a  function  of  en- 
ergy for  both  collimators.   The  two  inch  results  yield  a  smooth  curve.  As 
previously  mentioned,  the  half-inch  amplitudes  should  be  approximately  a 
factor  of  256  below  the  two-inch  amplitudes.   With  this  in  mind,  such  a 
curve  was  drawn  on  the  same  plot,  based  on  the  two-inch  data,  and  then 
the  half-inch  data  were  plotted.   The  scattering  of  the  half-inch  data 
about  the  theoretical  line  is  obvious,  the  worst  case  occurring  at  166 
keV,  where  the  discrepancy  is  nearly  a  factor  of  two.   It  is  believed 
that  these  discrepancies  were  caused  by  small  errors  in  positioning  the 
source  on  the  source  rack.   A  small  error  in  positioning  can  make  a  large 
difference  in  the  numerical  integration  of  the  point  source  data  for  the 
half-inch  collimator.   Figure  16  is  a  plot  of  the  photopeak  amplitudes  for 
the  point  Ce-139  source  as  a  function  of  the  distance  off  center.   It 
seems  that  the  maximum  of  the  distribution  occurs  somewhere  between  the 
assumed  center  and  1/4"  off  the  center.   The  numerical  integration  was 
done  from  zero  to  1.5  inches  and  gave  an  integrated  amplitude  of  0.000235, 
apparently  much  too  high  according  to  Figure  15.   A  second  integration  was 
done  using  1/4"  as  the  center  location,  this   gave  a  value  of  0.000147, 
a  factor  of  1.6  less  than  the  previous  calculation.   This  is  a  very  large 
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error  for  such  a  small  error  in  positioning.   Obviously  positioning  is 
very  critical  in  the  half-inch  experiments.   It  seems  as  i  f  these  types 
of  systematic  errors  are  present  in  at  least  three  of  the  half-inch  i 
perimcnts  and  hence  the  discrepancies  noted.   It  should  be  noted  that  the 
source  "dot"  was  approximately  0.125"  in  diameter  for  all  the  experiments, 
and  an  error  of  one  or  two  source  diameters  near  the  center  position  will 
effect  the  final  results  drastically. 

The  half-inch  collimator  response  functions  obtained  in  this  work 
are  presented  in  tabular  form  in  Appendix  D  together  with  their  appro- 
priate variances.   The  variances  obtained  in  the  half-inch  experiments  are 
not  as  good  as  those  obtained  with  the  two-inch  collimator,  and  usually 
run  around  10  percent.   However,  it  is  obvious  that  large  systematic  errors 
are  present  in  some  of  the  results,  and  in  at  least  one  case  makes  the 
statistical  errors  almost  insignificant.   The  results  are  presented,  how- 
ever, as  good  "first"  data  for  the  half-inch  collimator,  since  no  attempt 
had  been  made  previously  at  KSU  to  measure  the  half-inch  response. 

4.4  Conclusions  and  Suggestions 
The  integrated  point  source  method  is  capable  of  yielding  consistent 
results  which  are  both  accurate  and  precise  if  extreme  care  is  taken  in 
the  source  positioning.   This  is  supported  by  the  fact  that  the  two-inch 
collimator  results  obtained  in  this  study  agreed  favorably  with  previously 
obtained  data  which  were  accumulated  on  the  KSU  system  independently  of 
this  study  using  a  completely  different  plane  source  simulation  technique. 
The  half-inch  collimator  results  cannot  be  supported  directly  since  no 
previous  data  are  available.   However,  since  the  experimental  technique 
was  the  same  for  both  collimators,  it  is  believed  that  good  results  are 
obtainable  with  the  smaller  collimator  if  a  more  accurate  source  positioning 
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technique  is  devised. 

It  was  intended  to  accurately  calibrate  the  tv.7o-inch  collimator 
system  in  order  to  check  previous  experiments  and  provide  accurate  re- 
sponse functions  for  future  experiments;  this  goal  was  achieved.   It  was 
also  intended  to  obtain  first  data  on  the  half-inch  collimator  system  in 
order  to  observe  general  trends  and  relationships  between  the  two  systems. 
This  goal  was  also  achieved  since  some  basic  resolution  and  relative 
magnitude  trends  became  apparent  which  were  not  obvious  before  the  exper- 
iments were  initiated. 

A  major  disadvantage  of  the  integrated  point  source  method  is  the 
amount   of  data  which  must  be  accumulated  in  order  to  obtain  reasonable 
results.   In  this  study  approximately  750  spectra  were  obtained  over  an 
eight  month  period.   The  data  reduction  and  analysis  required  approxi- 
mately ten  hours  of  IBM  360/50  machine  time  and  many  more  hours  of  manual 
data  handling.   It  is  believed  that  comparable  results  are  obtainable 
using  other,  faster,  means  of  plane  source  simulation,  such  as  the  plane 
simulator  method  described  by  Baran  (1).   The  plane  simulator  method  would 
probably  be  better  for  the  half-inch  experiments  than  the  method  used, 
since  positioning  should  be  somewhat  easier,  although  not  trivial  by   any 
means . 

In  this  study,  better  statistical  results  could  have  been  obtained 
if  two  different  sets  of  sources  had  been  used.   The  sources  used  were  de- 
signed for  the  two-inch  collimator  experiments.   These  sources  were  not  in- 
tense enough  for  the  half-inch  experiments  because  of  the  factor  of  sixteen 
reduction  in  the  count  rates.   This  resulted  in  long  accumulation  times  for 
the  half-inch  experiments,  which  is  usually  undesirable.   If  the  sources 
used  had  been  sufficiently  intense  for  the  half-inch  collimator,  they  would 
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have  been  too  intense  foi-  the  two-inch  collimator  expel  Lmi  nls  and  would 
have  caused  large  dead  times  In  the  counting  system  which  could  have  in- 
validated the  data  because  of  rapid  spectral  shifts  associated  with  large 
dead  times. 

Quality  control  during  spectrometer  experiments  is  a  continuing 
requirement.   The  baseline  and  the  energy  per  channel  must  be  well  known 
during  the  experiments  or  the  accumulated  data  may  be  useless.   In  pre- 
vious studies  and  in  this  study,  so-called  quality  control  experiments 
were  performed  using  a  combination  Co-60-Cs-Ba-137  source  for  the  cali- 
brations.  These  experiments  were  done  approximately  every  two  hours  dur- 
ing an  experiment  day.   This  process  is  very  expensive  both  in  time  lost 
from  actual  experiments  and  in  computer  time  required  for  the  analysis  of 
the  quality  control  spectra.   It  is  suggested  for  future  work  that  a 
commercially  available  stabilizer  be  purchased  for  the  collimator-spec- 
trometer  system.   Such  a  feedback  device  could  maintain  the  energy  per 
channel  and  the  zero  intercept  at  preselected   levels.   This  would  greatly 
reduce  the  need  for  quality  control  experiments,  but  not  eliminate  the 
need.   In  this  author's  opinion,  the  savings  in  time  and  money  would  be 
significant. 

The  response  functions  presented  in  this  work  are  valid  only  for  the 
KSU  spectrometer  system  as  it  presently  stands.   It  is  known  that  as  the 
detector  ages   the  resolution  of  the  crystal  becomes  poorer  (1).   Event- 
ually the  response  functions  will  not  properly  represent  the  system.   Also, 
if  the  crystal  detector  is  replaced  for  some  reason  the  system  will  have 
changed  significantly  and  the  response  functions  presented  will  be  invalid. 
It  is  suggested  that  a  study  be  made  as  to  whether   or  not  existing  response 
functions  could  be  modified  to  reflect  these  system  changes.   The  modifi- 
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cations  of  course  would  have  to  be  based  on  experimental  data.  It  is 
believed,  however,  that  the  effort  required  to  accumulate  the  new  data, 
on  which  to  base  the  modifications,  would  be  less  than  that  required  to 
accumulate  a  new  set  of  response  functions. 
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7.1  APPENDIX  A 
Source  Strength  Calibration 

Three  methods  were  used  to  determine  the  source  strengths  of  the  nine 
isotopes  used  in  this  study.  Four  source  strengths  were  determined  from 
absolute  counting  experiments;  two  more  were  obtained  by  doing  comparative 
counting  experiments  using  National  Bureau  of  Standards  certified  sources. 
The  remaining  three  source  strengths  were  based  on  a  calibration  line  con- 
structed using  data  from  the  other  six  sources. 

Initially  it  was  hoped  that  all  of  the  sources  could  be  calibrated 
using  absolute  counting.   This  was  not  possible  since  the  sources  were 
originally  in  liquid  form;  consequently  if  all  of  the  liquid  were  dried 
on  the  planchet,  five  of  the  sources  would  have  been  too  "hot"  to  count 
in  the  absolute  counting  geometry  available.   No  reliable  method  for  dry- 
ing part  of  the  source  solution  on  the  planchet  and  keeping  an  accurate 
account  of  the  volumes  used  was  available  since  the  source  volumes  were 
typically  less  than  0.1  ml  and  buret  surface  tension  and  evaporation  be- 
came critical  in  the  volume  measurements.   Absolute  counting  was  therefore 
only  possible  with  the  four  sources,  Hg-203,  Sr-85,  In-Sn-113,  and  Ce-139. 
Mn-54  and  Cs-137  were  calibrated  by  comparing  the  photopeaks  of  the  exper- 
imental source  spectra  with  those  of  the  NBS  certified  source  spectra, 
both  pairs  of  which  were  accumulated  under  the  same  conditions.   Cd-109, 
Pb-210,  and  Zn-65  were  calibrated  using  a  technique  to  be  discussed  later. 

7.1.1  Absolute  Counting  Theory 

The  detection  efficiency,  T(E),  of  the  crystal  is  the  first  major 

parameter  of  importance.   For  well  defined  source-detector  geometries, 

this  efficiency  may  be  calculated  from  the  known  values  of  the  absorption 
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cross  sections,  x(E)   (20),  for  Nal.   Equations  have  been  derived  for 
several  cases  (7,20).   A  monoenergetic  point  source  of  radiation  located 
on  the  extended  axis  of  a  right  circular  cylindrical  detector  is  of  in- 
terest here.   Figure  17  shows  the  configuration  used  in  this  experiment. 


Point 
Source 


Detector 


Figure   17.        Source-detector  geometry   for   absolute  counting 
experiments . 
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Equation  (A-l)  is  the  analytical  expression  used  to  evaluate  T(E).   i 
integral  expressions  were  evaluated  numerically  using  a  96  point  Gaussian 
quadrature. 

rA 


T(E)  =  1/2  • 


[  1-exp  (-  T^„o)   ]   sinOdO  + 
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]  )  ]  sinOdO 
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A  =  tan  1    (  r^J—  ) 
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B  -  tan  1    (  ^   )  . 
ho 


The  calculated  absolute  detection  efficiency  obtained  from  the  above 
expression  is  the  probability  that  a  photon  of  energy  E  emitted  from  the 
source  will  interact  in  the  detector  with  a  finite  energy  loss.   Any  error 
to  be  expected  will  be  due  to  the  uncertainty  in  the  tabulated  values  of 
x(E)  used  in  the  calculation.   The  maximum  uncertainty  in  t (E)  is  at  the 
higher  energies,  greater  than  1  MeV,  but  it  is  still  less  than  10  percent, 
which  results  in  an  error  of  about  5  percent  in  T(E)  (7,20).   The  error  is 
reduced  considerably  at  the  lower  energies,  where  a  3"  x  3"  crystal  de- 
tector is  practically  opaque  to  gamma  radiation  of  energies  less  than  200 
keV. 

The  number  of  events  observed  in  the  pulse-height  distribution  is  re- 
lated to  the  emission  rate  of  the  source.   The  materials  surrounding  the 
detector  may,  however,  scatter  gamma  rays  into  the  crystal  which  tend  to 
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mask  this  relationship.   If  such  contributions  do  exist,  they  must  be 
removed  in  order  for  this  relationship  to  be  clear.   Since  these  contri- 
butions are  of  lower  energy  than  the  source  energy,  it  is  convenient  to 
perform  the  calculations  in  terms  of  the  photopeak  efficiency,    £  (E) . 
This  is  defined  as  the  probability  that  a  gamma  ray  of  energy  E,  emitted 
from  the  source  will  cause  a  count  in  the  photopeak  of  the  observed  spec- 
trum.  The  photopeak  represents  a  more  accurate  measure  of  the  emission 
rate  because  other  contributions  from  scattering  have  virtually  no  effect 
on  the  photopeak  in  comparison  to  effects  on  the  rest  of  the  spectrum. 

Direct  calculation  of  e  (E)  represents  a  formidable  problem  because 
of  the  large  number  of  multiple  processes  which  can  occur  in  the  crystal. 
Consequently,  it  is  convenient  to  use  the  following  expression: 

e  (E)  =  T(E)P(E)  (A-2) 

P 

where  T(E)  is  as  before,  and  P(E)  is  defined  as  the  fraction  of  the  total 
number  of  events  in  the  spectrum  which  appear  in  the  photopeak,  i.e.,  the 
peak-to- total  ratio.   The  peak-to-total  ratio  is  determined  experimentally 
by  careful  measurement  of  monoenergetic  sources  under  experimental  con- 
ditions which  reduce  the  scattering  effects  to  an  insignificant  level. 

The  emission  rate  of  a  monoenergetic  isotope  is  given  by  the  following 
expression  (7) : 

N 

N   =  2 (A-3) 

o     T(E)P(E)A(E)t  K        J 

where  N  is  the  number  of  gamma  rays  emitted  from  the  source  per  experiment. 
N   is  the  total  number  of  counts  under  the  photopeak,  and  A(E)  is  an  absorber 
correction  factor  used  to  account  for  the  attenuation  of  the  gamma  rays 
caused  by  the  insertion  of  a  beta  shield  between  the  source  and  the  detector  , 
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and  t  is  the  experiment  live  time. 

The  absorber  correction  factor  cannot  be  computed  directly  from  tab- 
ulated cross  section  data,  since  these  values  are  usually  reported  only 
for  highly  collimated  geometries,  it  must  therefore  be  measured  experi- 
mentally for  each  source  energy  used.   The  best  measure  of  A,  as  supported 
by  the  form  of  equation  (A-3 ) }is  the  ratio  of  the  total  number  of  counts 
in  shielded  and  unshielded  spectra. 

7.1.2  Counting  Procedure 

The  source-detector  assembly  was  placed  inside  a  large  counting  cav- 
ity with  inside  dimensions  of  32"  x  32"  x  32".   Figure  18  is  a  photograph 
of  the  partially  completed  cavity.   The  cavity  is  lead  lined  with  cadmium 
and  copper  grading  introduced  to  reduce  the  intensity  of  the  character- 
istic lead  X  rays  which  may  be  produced  in  the  lining.   This  type  of  shield- 
ing tends  to  minimize  the  scattering  effects,  which  is  desirable  in  any 
spectrum  analysis. 

The  source  was  placed  10  cm  above  the  crystal  center  by  allowing  the 
source  planchet  to  rest  on  a  plastic  mount  with  a  direct  opening  to  the 
crystal  face.   Thus,  a  standard  geometry  (7)  is  achieved.   The  detector  used 
was  a  Harshaw  Integral  Line  assembly  with  a  3"  x  3"  Nal(Tl)  scintillation 
crystal,  operated  at  +900  volts.   The  output  pulses  from  the  photomultiplier 
tube  are  fed  to  a  TMC  4096  multi-parameter  system. 

After  initial  planchet  loading,  the  source  was  counted  twice  under  the 
same  conditions,  except  that  one  spectrum  was  accumulated  with  the  beta 
shield.   Figure  19  is  a  photograph  of  the  detector  assembly  as  placed  in  the 
counting  cavity  with  the  beta  shield  in  place.   The  counting  procedure  was 
repeated  several  times  in  order  to  check  the  reproducibility  of  the  counting 
experiment.   After  the  source  runs  or  during  the  source  drying  procedure 
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background  runs  were  made.   These  runs  were  relatively  long  (120  rain.) 
so  as  to  insure  as  small  an  error  as  possible  in  the  background  deter- 
mination.  Again,  two  runs  were   completed   during   each   experiment 
day,  one  with  the  beta  shield  and  one  without  the  shield. 

7.1.3  Data  Analysis 

Figure  20  shows  a  typical  gamma-ray  spectrum  for  Hg-203.   This  par- 
ticular spectrum  was  obtained  without  the  beta  absorber  in  place  and  will 
be  used  to  indicate  how  a  typical  analysis  proceeded. 

Backscatter  and  other  wide  angle  scatter  contributions  are  much  re- 
duced in  the  absolute  counting  geometry  used,  but  are  still  visible  and 
must  be  removed  from  the  spectrum  since  they  are  not  part  of  the  response 
associated  with  the  source  energy.   The  broad  peak  in  channel  70  of  the 
Hg-203  spectrum  is  the  backscatter  peak.   The  characteristic  72  keV 
Tl  X-ray  peak  in  channel  40  and  its  associated  iodine  escape  peak  in 
channel  23  must  also  be  removed  since  they  also  are  not  part  of  the  crystal 
response  to  Hg-203  photons. 

The  Compton  distribution  on  either  side  of  the  backscatter  peak  is 
very  level  and  a  straight  line  may  be  drawn  between  them.   This  procedure 
is  illustrated  in  Figure  21.   The  procedure  is  a  little  more  difficult  in 
the  region  of  the  X-ray  and  the  escape  peak.   A  Compton  distribution  is 
still  visible,  but  it  is  not  constant  from  one  side  of  the  peaks  to  the 
other.   Also,  since  the  two  peaks  are  not  completely  resolved,  the  distance 
between  the  extremes  of  the  peaks  is  too  great  to  warrant  a  simple  averag- 
ing process.   Consequently,  since  the  Compton  distribution  is  slightly 
increasing  with  decreasing  energy,  a  least  squares  straight  line  was  deter- 
mined using  four  points  on  either  side  of  the  pair  of  peaks,  the  result  be- 
ing 
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Figure  21.   Illustration  of  spectrum  correction  technique,  Hg-203 
backscatter  peak. 


63 


where  i  is  the  channel  number  and  Y(x.)  is  the  calculated  count  rate  in 
channel  i. 

Having  subtracted  the  undesirable  contributions  from  the  "shielded" 
and  "unshielded"  spectra,  one  has  a  good  representation  of  the  corrected 
crystal  response  to  279  keV  gamma  rays  for  both  shielded  and  unshielded 
experiments.   The  two  types  of  spectra  are  corrected  in  the  same  manner. 
One  is  now  equipped  with  the  necessary  experimental  data  to  determine  the 
absolute  gamma-ray  emission  rate  of  the  source. 

The  first  step  is  to  fit  the  photopeak  to  a  Gaussian  distribution 
and  to  determine  the  area  under  the  photopeak.   As  stated  in  section  2.2, 
the  fit  is  accomplished  by  the  use  of  a  polynomial  least  squares  tech- 
nique.  The  photopeak  area  is  available  from  the  fit  parameters  as 

N   =  Y  (ttB  )1/2   .  (A-5) 

p      o   o 

The  total  area  is  determined  by  integrating  the  entire  spectrum,  i.e., 
totaling  the  counts  in  the  entire  spectrum.   The  peak-to-total  ratio 
is  given  simply  as 

P(E)   -  N  /N  (A- 6) 

P   t 

where  N  is  the  total  area.   The  beta  absorber  correction  factor  is  the 
only  undetermined  parameter.   Inspection  of  the  emission  rate  equation 
reveals  the  experimental  quantities  required  to  measure  A. 

A  =   (total  counts)  /(total  counts)  (A-7) 

where  s  refers  to  the  "shielded"  spectrum  and  u  refers  to  the  "unshielded" 
spectrum. 

The  standard  deviations  associated  with  each  experimental  quantity 
were  calculated  and  the  final  emission  rate  is  reported  with  a  percent 
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standard  deviation.   It  should  be  noted,  however,  that  the  standard 
deviations  reported  do  not  include  the  systematic  errors  associated  v.7ith 
the  experiments.   Evidence  is  presented  below  which  indicates  that  the 
systematic  errors  are  smaller  than  the  statistical  errors. 

In  order  to  check  the  validity  of  the  absolute  counting  technique 
described,  three  independent  experiments  were  performed  using  NBS  certi- 
fied sources.   The  counting  procedure  for  the  NBS  sources,  Mn-54  and 
Cs-137,  was  the  same  as  that  for  the  experimental  sources.   The  analysis 
of  the  data  was  done  in  the  same  manner  in  order  to  insure  consistency. 
The  NBS  Cs-137  source  was  counted  on  two  different  days  in  order  to  check 
the  reproducibility  of  the  experiments.   All  of  the  results  are  normal- 
ized to  the  same  count  rate  and  corrected  for  radioactive  decay.   The  re- 
sults are  shown  below. 

Table  I 
Verification  of  Absolute  Counting  Technique 
Source  Strengths,  (sec   ) 

Isotope NBS Measured 

Mn-54  1.75  x  104  +  1.0%         1.81  x  104  +  5.1% 

Cs-137         4.55  x  104  +  1.5%         4.48  x  104  +  3.5% 
Cs-137         4.55  x  104  +  1.5%         4.40  x  104  +  4.1% 


These  results  are  considered  to  be  adequate  proof  of  the  validity  of  the 
technique  used.  The  results  also  indicate  that  systematic  errors  in  the 
absolute  counting  experiments  were  small. 
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7.1.4  Comparative  Counting  Tccl tnlqi 
Comparative  counting  is  a  relatively  simple  technique.   Two  NBS 
certified  sources  were  available  for  use  in  this  experiment,  Mn-54  and 
Cs-Ba-137.   The  experiments  were  performed  in  collimated  geometry.   The 
source  planchet  was  placed  on  the  collimator  face  so  that  the  source  dot 
was  in  the  center  of  the  aperture.   Several  experiments  were  done  with 
both  the  NBS  sources  and  the  corresponding  experimental  calibration 
sources  under  the  same  conditions.   Each  spectral  photopeak  was  then 
fitted  and  the  photopeak  area  calculated  from  the  fit  parameters.   The 
photopeak  areas  were  then  averaged  for  each  source  and  the  appropriate 
standard  deviations  calculated.   The  emission  rate  for  the  calibration 
source,  was  then  calculated  from  the  NBS  source  strength  using  the  follow- 
ing expression. 


Ncal  _  ^BS 

^"*  "  "¥1  (A-8) 

where  N     is  the  certified  source  strength  for  the  given  isotope  and 

NBS   ° 
N     is  the  calculated  photopeak  area  for  the  NBS  source,  N     is  the 

calculated  photopeak  area  for  the  experimental  calibration  source,  and 

Ncal  is  the  required  emission  rate  for  the  calibration  source.   The  tech- 
o 

nique  is  valid  because  it  can  be  shown  that  the  photopeak  area  is  direct- 
ly proportional  to  the  emission  rate  (7,20). 

An  expression  for  the  standard  deviation  of  the  calculated 
source  strength,  derived  from  standard  propagation  of  error  techniques 
(2),  is  given  on  page  66: 
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NCal  NNBS 

cal    r   /  P  \         i     NBS  \2  _i  /   °    \2  ,  cal.?  , 

%        =  [   (  ^BS  }    (  V   }     (  >^  }   (°N   } 
o         N  °         N  p 

P  P 


NBS   cal 
_o p__.)2     NBS   2   1/2   _  (A_9) 


,.tnbsn2       n 

(N    )  P 


After  the  calculation  is  complete,  the  calculated  source  strength  for  the 
calibration  source  is  corrected  for  radioactive  decay. 

7.1.5   Calibration  Line  Technique 

SD(E)  is  considered  to  be  a  factor  which  accounts  for  collimator 
is. 

effects  in  the  experimental  spectra.   These  effects  are  primarily  mani- 
fested in  lip  penetration  and  wall  scattering  in  the  collimator  barrel. 

It  has  been  experimentally  shown  (1,11)  that  the  following 
equation  displays  an  approximately  linear  nature  with  energy. 


T(E) 


SR(E)  - 


K(l-e  u(E)t)   [  E1(b(E))  -  E1(b(E)secG)  ]   .    (A-10) 


T(E)  is  the  total  number  of  counts  in  the  normalized  spectrum  due 
to  an  isotope  emitting  photons  of  energy  E.    The  spectrum  is  normalized 
to  a  unit  source  strength  and  a  unit  live  time.   Mather  (11)  indicates 
that  S  (E)  is  unity  at  zero  energy.   K  is  a  proportionality  constant  which 
assures  this.   The  term   (  1-e      )  is  representative  of  the  crystal 
efficiency,  y (E)  is  the  total  attenuation  coefficient  for  gamma  rays  of 
energy  E,  and  t  is  the  crystal  thickness  (three  inches). 

The  angle   6   is  the  maximum  angle,  between  the  collimator  center 
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line  and  a  line  drawn  from  the  center  of  the  crysta]  face  to  the  oul 

edge  of  the  collimator  barrel.   E  (b(K))  is  the  first  order  exponential 
integral  for  b(E)  mean  free   paths  in  air  from  the  source  plane  to  the 
crystal  face  for  the  isotope  emitting  photons  of  energy  E.    The  term 
[  E  (b(E))  -  E  (b(E)secO)  ]  is  proportional  to  the  uncollided  flux  from 
a  disk  source  at  the  crystal  face. 

The  two  previously  discussed  source  calibration  techniques  pro- 
vide enough  information  for  the  calculation  of  six  points  on  the  S  (E) 
calibration  line.   These  six  points  are  then  fitted  with  a  straight  line 
in  the  least  squares  sense,  and  normalized  to  unity  at  zero  energy,  the 
result  being, 

S_(E)   =  6.719  x  10~4E  +  1.000 
R 

where  E  is  in  keV.   The  proportionality  constant  K  is  the  device  used  to 
normalize  the  calibration  line.   Figure  22  shows  the  resultant  calibra- 
tion line. 

With  the  calibration  line  at  hand,  the  procedure  used  to  calcu- 
late the  remaining  three  source  strengths  is  straight  forward.   S  (E)  is 
known  for  the  three  isotopes  under  consideration.   Equation  (A-10)  may 
then  be  solved  for  T(E),  the  total  number  of  counts  in  the  normalized 
spectrum,  i.e.,  the  ratio  of  the  total  count  rate  to  the  source  strength, 
Since  the  count  rate  is  known,  the  source  strength  may  be  determined 
directly. 

The  results  of  the  three  calibration  techniques  discussed  are 
shown  in  Table  II. 
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factor  as  a  function  of  energy  for  the  two-inch 
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Table   II 
Source  Strengths  of  Experimental  Isotopes 


Isotope 


E(keV) 


Technique 


Source  Strength 
(min-1) 


Pb-210 

47 

calib.* 

3.93 

X 

10 

+ 

5.6% 

Cd-109 

88 

calib. 

2.85 

X 

io6 

+ 

4.4% 

Ce-139 

166 

abs.ct.** 

3.44 

X 

io7 

+ 

1.1% 

Hg-203 

279 

abs.ct. 

2.33 

X 

io7 

+ 

5.8% 

In-Sn-113 

392 

abs.ct. 

2.35 

X 

io7 

+ 

3.0% 

Sr-85 

513 

abs.ct. 

4.30 

X 

io6 

+ 

2.5% 

Cs-137 

662 

comp.ct. t 

1.42 

X 

io8 

+ 

1.7% 

Mn-54 

840 

comp. ct. 

8.16 

X 

io7 

+ 

3.8% 

Zn-65 

1114 

calib. 

1.32 

X 

io8 

+ 

4.1% 

Calibration  line,  S, 


**  Absolute  Counting 
t   Comparative  Counting 


70 


7 . 2  APPENDIX  B 
Beta  Shield  Correction  Technique 

The  introduction  of  a  beta  shield  into  the  spectrometer  system  pro- 
duces two  primary  effects,  (1)  the  entire  gamma-ray  energy  spectrum  is 
reduced  by  some  constant  factor,  and  (2)  the  valley  between  the  photopeak 
and  the  Compton  continuum  tends  to  fill  up  because  of  small  angle  scatters 
in  the  beta  shield.   A  third  effect,  the  production  of  bremmstrahlung 
radiation  in  the  shield,  is  possible  but  is  considered  to  be  a  negligible 
effect. 

Source  energy  photons  which  reach  the  beta  shield  may  be  lost  from 
the  system  by  absorption  in  the  shield  or  by  backscattering  from  the 
shield.   It  is  necessary  to  derive  an  expression  for  the  backscattering 
cross  section;  this   is   done  by   integrating  the  Klein-Nishina 
differential  scattering  cross  section  over  the  back  hemisphere.   The  ab- 
sorption cross  section  is  then  added  to  the  backscattering  cross  section 
to  obtain  the  total  loss  cross  section.   The  derivation  is  outlined  below. 


2 

da  (6)  re 

— £ "  2- 

dfl     2  [  l+K(l-cos8)   ] 


2    K2(l-cos9)2 

1+cos  8  +  ,  ,W1 rr 

l+K(l-cos0) 


(B-l) 


2 
=  Klein-Nishina  Compton  scattering  cross  section,  (cm  /elec)/sr  (18) 

2 
where  K  =  E  /m  c   and  r   is  the  classical  electron  radius.   In  order  to 
o   o  e 

simplify  the  integration  process,  the  following  change  of  variables  is 
made. 

u  =  1  -  cos6 
du  =  sinG  d9  (B-2) 


7J 


This  leads  to: 


2     2  2  ? 

do  (u)     r     u  -2u+2     K  u 

"  du"       2   L  (1+Ku) 2     (1+Ku)3  J 


(B-3) 


KN 
inc,BS 


=  Klcin-Nishina  incoherent  backscattering  cross  section 

c     smGdO 
Tf/2   d6 


=  2tt 


da  (u) 
c 

du 


(B-4) 


du 


=  TTr 


2, 

u  du 

(1+Ku) 


_udu 
(1+Ku) 


2  +  2 


du     ,  2 
(Mu)2  +  K 


2, 

u  du 

(1+Ku) 3 


The  four  integral  expressions  in  equations  (B-4)  may  be  evaluated  separately 
and  then  recombined  to  yield  the  final  result.   Equation  (£-4)  is  rewritten  as 


Q™   mc  =  Trr2   (I  -2I0+2I„+K2I.) 
xnc,BS     e    12   3    4 


(B-5) 


where, 


rl-J 

'2 
1 

2, 

u  du 

(1+Ku) 2 

T2   = 

(2 
1 

udu 
(1+Ku) 2 

4 

'2 

1 

du 
(1+Ku)  2 

H' 

'2 

2, 
u  du 

i 

(1+Ku) 3 

(B-6) 
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With  aid  of  tables  of  integrals  (6)  ,  equations  $-6)  axe   evaluated  as 


I  =  *  +1  in  [(*!i)2  i  + 1 

Ll         K?  +  p  ln  L  k2K+l;   J  +  K2(K+1 


k2(K+1)(2K+1) 


T   _  _1_  ,   r  /2K+1.  2  i L_ 

X2  "  2K2       K+l  ;   J    K(K+1)(2K+1) 


X3    (K+1)(2K+1 


•*■* 


i  /2K+1. 

ln(R«r) 


2K 


+ 


K(3K+2) 


(2K+1XK+1)    2(2K+1)2(K+1)2 


(B-7) 


Combining  the  four  integrals  according  to  equation  (B-5) ,  we  write  the 
expression  for  the  backscattering  cross  section. 


KN 
a..        „(K) 


inc,BS 


=  irr  S 


1   r  4K2+5K+2 
Kz 


(K+l) (2K+1) 


2K+1 
2 (K+l)  ,   ,2K+1N1  ,  ln^K+1  ) 
K~  ln  (^+l)]  +  K 


8K  +9K+2 


2(2K+1)2(K+1)2 


cm  /elec 


(B-8) 


The  shield  was  made  of  lucite  (polymethyl  methacrylate,  C^H^O  ).   For 


computational  ease,  the  units  should  be  converted  from  cm  /electron  to 

cm  /g 
cm2/elec 


0  O  /i       f) 

cm  /gm,  using  the  conversion  factor  f  =  0.32501  x  10    cm  /g    ,  where 
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f  = 


N   Z 
a 


A  -   100.06256  graras/g-mole 

Z   =  54   electrons/molecule 

24 

N  =  0.60225  x  10   molecules/g  -mole. 

3 

For  the  five  isotopes  where  the  corrections  are  necessary,  lucite  displays 
an  appreciable  photoelectric  absorption  cross  section  only  for  the 
lead-210  energy.   At  this  energy   the  cross  section  was  evaluated  as  (18) 

o-Lucite  (Pb-210)  =  0.0185  cm2/g  • 

T,t 

For  mercury-203   the  photoelectric  cross  section  would  increase  the  total 
loss  cross  section  by  less  than  0.25  percent.   Above  the  mercury-203  en- 
ergy, no  photoelectric  cross  sections  are  available  for  lucite (18).   The 
total  loss  cross  sections  are  as  follows. 

Table  III 
Total  Beta  Shield  Loss  Cross  Sections 


Isotope 


Energy  (MeV) 


Cross  Section  (cm  /g) 


Pb-210 

Hg-203 

In-Sn-113 

Cs-137 

Zn-65 


0.047 
0.279 
0.392 
0.662 
1.114 


0.1023 
0.0411 
0.0337 
0.0243 

0.0170 


The  cross  sections  are  converted  to  gamma-ray  attenuation 
coefficients,  u  (cm  ),  by  multiplying  them  by  the  shield  density  which 
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3 
was  measured  and  found  to  be  1.17  gm/cm  .   Multiplying  again  by  the  shield 

thickness,  1/4",  yields  the  number  of  mean  free  paths,  yx.   The  whole  spec- 

yx 
trum  is  then  multiplied  by  e   to   account   for  the  shield  reduction 

effect. 

The  filling  of  the  valley  between  the  photopeak  and  the  Compton  tail 
is  assumed  to  be  a  minor  effect  (see  section  4.2).   It  is  believed  that 
a  technique  for  correcting  this  filling  effect  could  be  developed,  but  such 
a  development  is  beyond  the  scope  of  this  work. 

In  order  to  check  the  accuracy  of  the  multiplicative  correction 
applied,  data  accumulated  on  the  TMC  4096  multiparameter  system  were  used. 
The  data  consisted  of  both  "shielded"  and  "unshielded"  spectra  for  the  beta 
emitting  isotopes.   The  table  below  compares  the  total  area  of  the  "un- 
shielded" spectra  to  those  of  the  corrected  "shielded"  spectra. 


Table  IV 
Comparison  of  Beta  Shield  Corrected  Spectra 


Isotope     Energy  (MeV)   Cross  section  (cm  /g)    e     Percent  Deviati 


on- 


Pb-210        0.047            0.1023  1.0997  -5.3 

Hg-203        0.279            0.0411  1.0389  +2.0 

In-Sn-113      0.392            0.0337  1.0318  -0.4 

Cs-137        0.662            0.0243  1.0228  -0.2 

Zn-65         1.114            0.0170  1.0159  +1.2 

*  (corrected  -  unshielded)  /unshielded 

The  correction  seems  to  be  very  good  in  the  middle  energy  range,  poor 
in  the  lower  energy  range  and  adequate  In  the  higher  energy  range.   The 
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apparent  loss  of  area  in  the  lower  energy  i  inge  could  be      d  by  tl 

absorption  of  small  angle  scattered  photons  in  the  crystal  containment 
can.   It  is  also  indicated  (18)  that  the  uncertainties  in  the  scattering 
cross  sections  may  be  as  high  as  3  percent,  and  the  uncertainties  in  the 
absorption  cross  sections  may  be  as  high  as  10  percent.   With  these 
possible  errors  in  mind,  it  is  believed  that  the  beta  shield  correction 
is  adequate. 
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7.3  APPENDIX  C 
Spectrum  Fitting  Technique 

The  most  important  part  of  the  gamma-ray  spectrum  is  the  photopeak. 
The  photopeak  amplitude  is  proportional  to  the  source  strength  of  the 
isotope,  as  is  the  photopeak  area.   The  photopeak  is  less  susceptible  to 
spectral  abberations  such  as  small  angle  scatters  and  bremmstrahlung; 
it  is  therefore  the  "cleanest"  part  of  the  spectrum.   Consequently   the 
photopeak  is  capable  of  yielding  information  about  the  gamma-ray  emitting 
isotope  under  study. 

As  previously  mentioned,  the  photopeak  is  approximately  Gaussian  in 
shape.   The  Gaussian  shape  has  a  physical  basis  in  that  the  interactions 
in  the  crystal  are  statistical  in  nature  and  the  detector  and  PM  tube 
operate  in  a  statistical  manner.   The  photopeak  is  described  then  as, 

Y(x.)=  Yo  exp["(xi"Xo)   ]    .  (C-l) 

B 
o 

where  Y(x.)  is  the  count  rate  in  channel  x..   Y  is  the  maximum  photopeak 

1  X      O  r      r 

amplitude,  located  at  pulse  height  x  .   B  is  the  photopeak  width  par- 
ameter equal  to  twice  the  variance  of  a  normal  Gaussian  distribution.   The 

fitting  parameters  to  be  solved  for  then,  are  Y  ,  x  and  B  . 
°   r  o   o      o 

Equation  (C-l)  is  non-linear  and  the  least  squares  fit  of  the  data 
may  be  accomplished  in  a  non-linear  manner  by  using  an  iterative  process. 
This  type  of  process  is  not  very  economical,  however,  if  one  has  a  large 
amount  of  data  to  be  analyzed  on  a  computer.   An  alternative  is  to  linear- 
ize the  problem  by  transforming  the  equation  into  the  log  domain,  which 
results  in  a  quadratic  equation  that  may  be  fitted  by  using  standard  poly- 
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nomial  least  squares  techniques.   Comparisons  between  the  two  meth 
were  made  and  the  results  were  essentially  the  same,  with  the  Latter 
method  requiring  much  less  computer  time  than  that  required  by  tin 
iterative  method. 

The  basic  input  to  the  code  used  for  the  spectrum  analysis  is  the 
spectrum  itself,  the  number  of  channels  of  data  (400),  and  the  fitting 
limits.   The  fitting  limits  are  very  important,  since  these  two  parameters, 
the  upper  and  lower  limits,  define  between  what  channels  the  photopeak 
lies  and  therefore  where  the  problem  is  to  be  solved. 

It  has  been  standard  procedure  (7)  to  base  the  limits  on  an  es- 
timate of  the  photopeak  height.   The  lower  limit  is  defined  as  that  chan- 
nel in  which  Y   is  reduced  by  approximately  one  third,  the  upper  limit  is 
that  channel  in  which  Y   is  reduced  by  approximately  one  half.   The  limits 
are  chosen  in  this  manner  because  small  angle  scatters  and  crystal  non- 
linearities  cause  the  lower  part  of  the  photopeak  to  deviate  from  the 
Gaussian  shape.   The  limits  are  not  'symmetric  because  small  angle  scatters 
affect   only  the  lower  energy  side  of  the  photopeak.   This  method  of 
choosing  the  limits  results  typically  in  a  fifteen  channel  fit  of  the 
photopeak  using  nine  channels  above  x  and  six  channels  below. 

7.3.1  Outline  of  Problem  Formulation  and  Solution 
Let  L  represent  the  lower  fitting  limit  and  U  represent  the  upper 
limit.   The  least  squares  problem  is  represented  in  the  log  domain  as 

AS   =  B 
where  A  is  the  coefficient  matrix,  S  is  the  solution  vector  and  B  is  the 
constants  vector.   A  is  a  symmetric  matrix  and  is  formed  as  shown  on 
page  78. 
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All  -  J  (1-L)"  Wi 

A12=J   (i-L)3W.  -A21 
i=L 

A^  =  J      (i-L)2  W.  =  A   =  A22 

U 

A23  =  JL  (i-1}   \  =  A32 

U 

A33  =  I     Wi  *  (C~3) 

i=L 

where  W.  are  the  appropriate  weighting  factors,  as  previously  discussed 
in  section  2.2,  in  the  log  domain.  The  constants  vector  is  set  up  in  a 
similar  manner,  as  shown  below. 


v      3-1 
B  =  I    (i-L)J  3   z.W   ;  j  =  1,2,3.  (C-4) 

J   i=L 


where  z.  is  the  natural  log  of  Y(x.)  and  W.  is  again  the  weighting  factor, 

The  solution  vector  is  obtained  by  premultiplying  both  sides  of 
equation  (C-2)  by  the  inverse  of  the  coefficient  matrix, 

S  =  A_1B  (C-5) 


Sl  =  "1/Bo 

S  =  2x  B 

2  o  o 

S  =  In  Y  -(x  /B  )2   .  (C-6) 

3  o   o  o 
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The  fitting  parameters  are  then  given  as 


B   =  -1/S 
o         i' 

S2 
x   =   -1/2  -f      , 

Sl 

x  2 

Y   =   exp  (-?-  +  S„).  (C-7) 

O  D  J 

o 


The  coefficient  matrix  is  inverted  by  using  a  bordering  method  which  is  both 
fast  and  accurate  for  the  3x3  matrices  under  consideration. 

If  F.  represents  the  calculated  photopeak  height  in  channel  i,  the 
standard  deviation  of  the  fit  is  given  as 

r  Xft*-'i>2  ,1/2  (C"8) 

a  =    [  i-k J 

(U-L+l) 

The  extremes  of  each  fit  are  checked  in  order  to  insure  a  good  fit.   If 
(Y  -  F_ )  is  greater  than  1.5  a,  L  is  replaced  by  L+l  and  the  entire  fit 
is  redone.   If  (Y  -  F  )  is  greater  than  1.0a,  U  is  replaced  by  U-l  and 
the  problem  solved  again.   This  checking  procedure  was  designed  by  Heath 
(7)  and  is  considered  to  be  an  adequate  criterion  for  obtaining  a  good 
fit. 

The  photopeak  area  is  obtained  by  integrating  the  Gaussian  equation 
from  -°°  to  +  <=°  as  shown  below. 


N  = 
P 


Y(x)dx  =  Y  (ttB  )1/2   .  (C-9) 

o   o 


The  peak-to-total  ratio  is  then  immediately  available  as 


p(E)   =  N  /N    .  (C-10) 

P   t 
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The  matrix  formulation  of  the  least  squares  problem  provides  a 

straightforward  method  for  obtaining  the  variances  of  the  fit  parameters. 

The  variances  of  S.  are  taken  from  the  diagonal  elements  of  A   ,  i.e., 

A..   is  the  variance  of  S..   The  variances  of  the  fit  parameters  are  then 
n  1  r 

immediately  available  by  applying  standard  propagation  of  error  techniques 
(2). 

For  more  detailed  information  on  the  fitting  techniques  discussed, 
the  reader  is  referred  to  Heath   (7,8). 
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7.4  AITEND]X  D 


Tabulated  Response  Functions 

The  response  functions  presented  here  are  all  normalized  to  a  unit 
uncollided  flux  per  steradian,  such  that  the  units  are 

/    2        —1 
(counts/sec)   /  (cm  sec  sr) 

All  the  responses  and  their  associated  variances  are  tabulated  on  an  en- 
ergy scale  of  3.5  keV/channel,  with  zero  energy  occurring  at  the  beginning 
of  channel  two.   Channel  one  is  not  part  of  the  response  function  or  its 
variance.   The  variances  reflect  all  the  mathematical  operations  perform- 
ed on  the  response  functions  and  are  normalized  in  the  same  manner.   Each 
listing  is  titled  by  the  energy  of  the  response,  the  collimator  aperture 
used  and  whether  or  not  the  data  are  response  functions  or  their  assoc- 
iated variances. 

Each  page  contains  one  complete  response  function  or  response 
function  variance,  i.e.,  400  channels  of  data.   Two  hundred  channels 
appear  on  the  left  half  of  the  page  under  the  title  and  are  read  from 
left  to  right.   The  remaining  two  hundred  channels  of  data  appear  on  the 
right  half  of  the  page,  again,  read  from  left  to  right. 
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ABSTRACT 

Using  a  collimated  sodium   iodide  scintillator  wi  th  assoc- 
iated electronics  and  a  multichannel  pulse-height  analyzer,  pulse- 
height  spectra  were  measured  for  the  collimator  viewing  monoenergetic 
plane  sources  of  gamma  radiation.   These  "response  functions"  (ratios 
of  the  counts  per  pulse-height  unit  per  unit  time  to  differential  gamma- 
ray  flux  per  steradian  at  the  source  plane)  were  measured  for  nine 
gamma-ray  energies  between  46  keV  and  1.114  MeV,  and  for  two  collimator 
aperture  sizes  (two-inch   diameter  and  one-half  inch  diameter).   A 
point  source  superposition  technique  was  used  to  simulate  the  infinite 
plane  source. 

The  two-inch  collimator  response  functions  were  compared  with 
previously  measured  responses  on  the  same  system.   The  two  measurements 
were  made  independently  of  each  other  and  used  two  different  techniques 
for  simulating  the  infinite  plane  source.   The  response  functions  agreed 
favorably.   No  such  direct  comparison  was  possible  with  the  half-inch 
collimator  response  functions  since  no  previous  data  were  available. 

The  source  strengths  of  the  point  sources  were  determined  using 
three  different  methods.   An  absolute  counting  technique  was  used   for 
four  of  the  calibrations;  comparative  counting  was  used   for   two  of  the 
determinations.   National  Bureau  of  Standards  certified  sources  were  used 
for  the  comparative  method.   Three  of  the  determinations  were  based  on 
intercomparisons   of  data  with  data  from  the  other  six  measurements. 


